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Abstract 

The inherent properties of specific physical systems can be used as a metaphore for investigating 
the behavior of complex networks. This insight has already been put into practice by previous 
works, e.g., studying the network evolution in terms of phase transitions of quantum gases or 
representing distances among nodes as they were particle energies. In this work, we show that 
the emergence of a scale-free structure in complex networks can be represented in terms of a 
quantum-classical transition for quantum gases. In particular, we propose a model of fermionic 
networks that allows to investigate the network evolution and its dependence with the system 
temperature. Simulations, performed in accordance with the cited model, clearly highlight the 
separation between simple random vs. scale- free networks, with full correspondence with classical 
vs. quantum regions for quantum gases. 

PACS numbers: 89.75.-k, 89.75.Hc, 05.65.+b 



marco.javarone@diee.unica.it 


armano@diee.unica.it 





1 



I. INTRODUCTION 



Many natural and man-made complex systems can be modeled by complex networks, 
i.e., biological cells and their interactions, social dynamics, world wide web, and every other 
system containing a notable number of interacting elements (Ref. [I] [2| [2] ) • In general, com- 
plex networks allow to represent systems at a non-equilibrium state, meaning that new nodes 
and/or new links can be added over time during an evolutionary, often irreversible, process. 
In Ref. PP, Barabasi illustrated the central role of statistical mechanics for analysing these 
networks. One of the first models for random graphs dynamics has been developed by Erdos 
and Renyi, see Ref. [?]. This kind of networks, called ER graphs, have low clustering coef- 
ficient and a binomial degree distribution of nodes, converging to a Poissonian distribution 
for large number of nodes. Later on, in their WS model, Watts and Strogatz Ref. [5] in- 
terpolated ER graphs with a regular ring lattice, achieving networks characterized by small 
average shortest path length, high clustering coefficient, and homogeneous degree distribu- 
tion of nodes. Observing that in many real networks the degree distribution of nodes follows 
a power-law, Barabasi developed the BA model Ref. [T], defining the concept of scale-free 
networks. The power-law equation is the following: 

P(k) ~ kT 1 (1) 

where k is the node degree and 7 is a parameter usually in the range (2,3). A notable 
consequence is that, in these networks, only few nodes (called hubs) have a big number 
of links. The evolutionary process of networks has been studied by many authors, often 
resorting to theoretical physics. For the sake of brevity, let us cite only few. Bianconi 
and Barabasi Ref. [6] compared Bose-Einstein Condensation phenomena (BEC) to winner- 
takes-all policies; Bianconi Ref. [7] discussed the symmetric construction of bosonic and 
fermionic networks, asserting that the former networks are scale-free and the latter are 
growing Calay trees. Kriukov et al. Ref. [H] developed a geometric framework to study the 
structure and functions of complex networks, interpreting edges as non-interacting fermions 
whose energies are hyperbolic distances between nodes. Shen Yi et al. Ref. [9] discussed an 
inverse approach to network evolution defining a relation with their model and Fermi-Dirac 
statistics. Baronchelli et al. Ref. [10] defined a framework using bosonic reaction-diffusion 
processes, with the aim of analysing dynamical systems on complex networks. Perseguers et 
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al. Ref. [TT] developed a model of quantum complex networks, drawing a link between the 
field of complex networks and that of quantum computing. 

In this work we propose a theoretical model of network evolution based on Fermi-Dirac 
statistics, obtained by mapping complex networks to fermionic gases. The proposed model, 
which can be seen as dual with respect to bosonic networks (Ref. [0]), shows that the emer- 
gence of a scale-free structure can be represented in terms of a quantum-classical transition 
for quantum gases. In particular, at low temperatures networks are scale-free, whereas at 
high temperatures they become simple random networks. As a result of several simula- 
tions, we computed a separation line between scale-free and simple random regions, in a 
temperature-density diagram. Notably, this line performs a separation that mimics the one 
observable between classical and quantum regions of quantum gases. The remainder of the 
paper is organized as follows: Section [IT] gives a brief introduction to quantum statistics. 



Section III describes some existing models of networks as non-interacting particle systems. 



Section [TV] introduces the proposed model of network dynamics, based on Fermi-Dirac statis- 
tics, and shows the results of the corresponding simulations. Conclusions (i.e., Section |V| 
end the paper. 



II. QUANTUM STATISTICS 

Statistical mechanics assumes a central role when dealing with systems composed by many 
particles, the underlying assumption being that particles are identical and indistinguishable. 
Moreover, their quantum energy levels are extremely closely spaced, with a cardinality much 
greater then the number of particles. Energy levels can be grouped in bundles with the ap- 
proximation that levels in the same bundle have the same energy. Particles with a symmetric 
wave function, called bosons, obey Bose-Einstein statistics, whereas particles with an an- 
tisymmetric wave function, called fermions, obey to Fermi-Dirac statistics -see Ref. [T2] . 
Given a system with N particles of the same type, we can build an N-body wave function 
having in an a state a number of particles n a , called occupation number, such that: 



0,1,.. .,00 bosons 
n a = { (2) 

0, 1 fermions 



and J2 a n a = N. Considering a gas composed by N bosons, the number of microstates in 
such macrostate is computable by the following equation: 

n t = n,(^i±M (3) 

mlgtl 

with gi representing the i-th. bundle. The distribution of particles follows the Bose-Einstein 
statistics: 

< = -^f— (4) 

e k b T _l 

where denotes the energy of the i-th bundle, \i the chemical potential, and fc& the Boltz- 
mann constant. In the event that a gas is composed by fermions, we must consider also the 
Pauli exclusion principle. Overall, the number of microstates is computed by the equation: 

The distribution of particles follows the Fermi-Dirac statistics: 

4 = -Jt~ («) 

e k b T +1 

Both these distributions approach the classical behaviour in proximity of the high-temperature 
limit, showing a quantum-classical transition. Such phenomenon occurs when particles oc- 
cupy excited states sparsely. In particular, considering the thermal wavelength A and the 
density p, the following conditions can be stated: 

pA 3 ^> 1 classical regime 

(7) 

pA 3 1 onset of quantum effects 

The classical regime is described by the Maxwell-Boltzmann distribution. In particular, 

e j_ 

with Z = j 9j e kbT partition function, we write: 

nf = ^9ie-& (8) 



III. NETWORKS AS PARTICLE SYSTEMS 



Let us give a synthetic overview of some existing models of network dynamics based on 
quantum statistics: 
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A. Bosonic Networks 



In this model, Bianconi and Barabasi Ref. [6] compared network evolution to phase 
transition of bosonic gases. Two main structures, i.e., fit-get-rich and winner-takes- all, are 
identified as two different phases at low temperatures. In this model, each node is interpreted 
as an energy level and each link as a pair of particles. Starting from a fitness parameter 77, 
energy is computed according to the following equation: 

e = --logT7 (9) 

with /3 = Here, the fitness parameter r\ describes the ability of each node to compete for 
new links. In particular, for the z-th node, the probability of connection with new nodes is 
proportional to: 

(10) 

with hi degree of the z-th node. Notably, new nodes tend to link with pre-existing nodes 
having high values of (77, k). The generation of a scale-free network in the fit-get-rich phase 
is characterized by Equation ([I]) and entails the presence of few nodes with a high degree 
connected to many others with low degree. In a bosonic gas, when the temperature decreases, 
particles aim to occupy lower energy levels until T ~ 0, then Bose-Einstein condensation 
takes place. In this model, as T decreases, many particles move to lower levels while keeping 
the corresponding particles at upper levels. In doing so, links concentrate on few nodes, until 
they condensate in the winner-takes- all phase, characterized by the fact that only one node 
predominates. In Ref. [7j, Bianconi discussed the differences between bosonic and fermionic 
networks, showing that the former are scale-free, whereas the latter can be represented by 
Cayley trees. 



B. Fermi-Dirac Statistics of Complex Networks 

In this work Ref. [9], the authors proposed an inverse approach to network evolution, 
starting from a random network with an average degree equal to 2m. Using an illness 
model, at each step, one randomly-selected node becomes "ill", then its illness causes a 
loss of links. An illness parameter I is defined and assigned randomly to each node. The 
probability for each node to lose a link depends on / and on the node degree. Hence, the 
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number of links decreases with time. Each node has an energy defined by: 



-log/ (11) 



with (3 = Y?p. A link between two nodes corresponds to a couple of non-interacting particles 



placed in two energy levels whose value is computed by Equation 11 When a new node 
joins the network, a new energy level is added, with m new particles, and other m links 
are deleted, due to the illness. Authors showed that the behavior of such system can be 
approximated by a Fermi distribution, so that the network can be seen as a Fermi gas. 



IV. FERMIONIC NETWORKS 



Let us now introduce a novel proposal for modelling network dynamics, based on Fermi- 
Dirac statistics. Given a network with n nodes and v edges, say E = (n, v), let us represent 
each link as a particle and each node as a degenerate bundle of energy levels. Each energy 
level corresponds to a potential-link in the network, which is drawn only when at least one 
particle exists on it. The maximum number of links in a network is L = n(n — l)/2. Hence, 
in the limit case of a fully-connected network, the number of links is v — L. The number 
Qi of available states in the 2-th bundle is much bigger than the number p^ of particles, i.e., 
9i ^ Pi- The i-th level has an energy e i5 which corresponds to that of the hosting bundle. 
This value can be assigned randomly or depends on a property of the system, e.g., a fitness 
parameter r\ or any other function deemed relevant (with the trivial constraint that it must 
be computable for each node of the network). In the proposed model, lower bundles contain 
more energy levels. In particular, the first bundle has n — 1 levels, the second has n — 2 levels, 
and so on and so forth. Note that the link Zy, which connects nodes i and j, is represented 
only by a single energy level, i.e., ey, which in turn is contained in the i-th bundle (under 
the assumption that the i-th bundle is deeper than the j-th). In doing so, the last node "y" 
has no potential-links, although it can be linked in the event that a particle stays at the e xy 
level, with x corresponding to one of the other nodes. Figure [T] shows an example aimed at 
graphically illustrating these concepts. 
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FIG. 1. On the left, a small network with 6 nodes and 5 links (each node is labeled by a 
letter). On the right, the corresponding fermionic model: a bundle gt is assigned to each node, 
each potential-link has its proper energy level, and real links are represented by particles. 

A. Modelling Network Evolution 

Let us consider a dynamical network E = (n(t),v(t)), with n(t) of nodes and v(t) edges, 
both time- varying (e.g., the world wide web or specific social networks). According to the 
Fitness Model Ref. j6], when a new node with fitness r\ joins the network, pre-existing nodes 
compete to connect with it. 

For every node, a bundle is defined -whose energy is computed with Equation ([9]). The 
relative position of each bundle depends on the value of its energy, so that deeper bundles 
incorporate more available states. For each bundle, the number of available states increases 
if a new node with higher energy is added. Hence, nodes represented as bundles with 
low energy have more probability to gain new links. With t generic point in time, let the 
connectivity of the i-th node added at time ti be fcj(ej, i, Every node generates m links 
when added to the network, hence m new particles are assigned to lower bundles. We also 
assume that, as time passes and other nodes join the network, the connectivity of the i-th 
node increases, in agreement with the scale-free model Ref. [13] . according to the following 
power law: 

k l (e l ,t,t l )=m(^y^ (12) 
where /(e^) is an energy dependent dynamic exponent, defined as: 

f(e) = e"* (13) 

Considering the ability of particles to jump between energy levels, at high temperatures, the 
particle distribution follows the classical Maxwell-Boltzmann and particles are distributed 
among the available states -see Equation d8l). As temperature decreases, many particles 
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move to lower energy levels so that when new particles join the system, most likely they will 
be located at low energy levels too -see Figure |2j 



Simple Random Network 




FIG. 2. On the left, from top to bottom, the evolution of a network with 10 nodes and 9 links 
from a simple random network to a scale-free network. On the right, their corresponding fermionic 
models, which result from a cooling process that pushes particles to low energy levels, contained 
in two bundles. 

During this process, few nodes gain new links. Hence, when new nodes join the network, 
they have a higher probability to be selected, as their hi is increased. Given the number 
of particles, it is possible to compute the Fermi Energy Ef of the system as the energy of 
the bundle containing the last particle at T — > 0. Hence, as temperature decreases, and 
assuming that the number of particles approximates the number of bundles, the winner- 
takes-all phase takes place. On the other hand, in the event that the number of particles 
grows, a fit-get-rich phase emerges (see also Ref. [6]). We mapped the temperature to the 
scale parameter of the gamma distribution used to perform a weighted random selection of 
nodes. This selection is executed when a new node joins the network and generates m new 
links with pre-existing nodes. According to the fermionic model previously described, it is 
a weighted random selection, as nodes represented as bundles with more available states 
have more chances to win. Finally, as shown in the following, it is possible to draw a rough 
separation in a T — p plane between the simple random and the scale-free region. 
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B. Simulations of Fermionic Networks 



The proposed fermionic model has been tested with many simulations. In particular, we 
generated networks of different size, varying the value of m and the system temperature T. 
Figure [3] shows three networks, with n = 50 and m = 2, generated at different temperatures. 
At T = 0.5 and T = 3 the networks are scale-free, whereas at T = 10 a simple random 
network is generated. Their degree distributions are represented in the P(k) diagram. For 




FIG. 3. The figure shows three fermionic networks, generated at different temperatures, as well 
as their degree distributions P(k), with n = 50 and m = 2. a) A fermionic network generated at 
T = 0.5 (coloured in green). As shown in the diagram, its P(k) follows a power-law approximated 
by 7 = 3. b) A fermionic network generated at T = 3 (coloured in blue). As shown in the diagram, 
its P(k) follows a power-law approximated by 7 G (2,3). c) A fermionic network generated at 
T = 10 (coloured in magenta). It is not scale-free, as its P(k) does not follow a power-law. For the 
sake of reference, two referential curves for scale-free networks are also shown in the P(k) diagram: 
7 = 2 and 7 = 3, see Equation ([l]), are represented by a black and red dotted line, respectively. 

each simulation an approximated temperature of transition has been computed, from a scale- 
free structure to a simple random network. Figure [4] shows some results achieved using a 

9 



number of nodes from 250 to 2000, with different values of m. The temperature of transition, 




FIG. 4. Results of the degree distribution P(k) for some fermionic networks generated during the 
simulation, a) Fermionic network with n = 600 and m = 4. At temperatures T > 40 this network 
loses its scale-free structure, b) Fermionic network with n = 1200 and m = 2. At temperatures 
T > 70 this network loses its scale-free structure, c) Fermionic network with n = 1500 and m = 3. 
At temperatures T > 75 this network loses its scale-free structure, d) Fermionic network with 
n = 2000 and m = 2. At temperatures T > 110 this network loses its scale-free structure. Two 
referential curves for scale- free networks: 7 = 2 and 7 = 3, see Equation Q, are represented by 
black and red dotted lines, respectively. 



from a scale-free to a simple random structure, increases with the number of nodes. It is 
worth noting that these temperatures can be interpolated, in the T — p plane, by the curve 
pA 3 , with p number of particles (i.e., links of the network) and with A = aT~^ wavelength. 
The value a = 0.2085 best interpolates the temperatures - see Figure [5j This curve, for 
values less than the scale factor a, mimics the curve that draws a rough separation between 
classical and quantum regime for quantum gases -see also Equation (|7j). 
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FIG. 5. Simple Random Networks and Scale- free Networks regions in the T — p plane. The dotted 
line is pA 3 , with p density and A thermal wavelength. This line represents a rough separation 
between the two regions. The former corresponds to the classical region and the latter to the 
quantum region for quantum gases. Green O represent the transition values (T, p) (with an error 
of ±5) of the simulated fermionic networks from scale-free to simple random structures. 



C. Discussion 



Fermionic networks show that the emergence of a scale-free structure can be represented 
as a quantum-classical transition for quantum gases. In particular, a scale-free network 
correspond to a fermionic gas approximated by the quantum regime at low temperatures. 
On the other hand, a simple random network corresponds to the same gas in classical regime 
at high temperatures. Similar considerations about the connection between classical random 
and scale-free networks have been proposed in Ref. [8]. The authors show that, in the cold 
regime, their network is scale-free, but as the temperature increases, the network loses 
its metric structure and its hierarchical heterogeneous organization -becoming a classical 
random network. Considering that many real complex networks are scale-free and others 
have not this structure, see Ref. [H], we deem that the proposed fermionic model can be 
considered a good candidate for representing their evolution, at low and high temperatures, 
respectively. 
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V. CONCLUSIONS 



In this work, we represent complex networks as quantum gases, thus defining a fermionic 
network model. Using this model, we show that network evolution is a temperature- 
dependent process. We found that it is possible to draw a separation between simple random 
region and scale-free region, in the T — p plane, with a curve similar to that used to sepa- 
rate classical and quantum regions for quantum gases. The transition from simple random 
to scale-free networks mimics a cooling process of a physical system achieving equilibrium, 
despite the non-equilibrium nature of network evolution. During this process, the transition 
from classical to quantum regime (or vice-versa for heating processes) takes place. Finally, 
we believe that many real complex networks can be mapped to quantum systems, at low 
or high temperatures depending on their structure. As for future work, we are planning to 
develop a model in which the temperatures can vary during network evolution. 
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